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ABSTRACT 

In  the  present  paper  we  extend  our  recant  work  on  the  continuous 
single  nodule  design  problea  to  the  multiple  nodule  cue.  It  is  aseuaed 
that  there  is  a  fixed  cost  associated  with  each  additional  nodule  used 
to  solve  the  problem,  The  Kuhn- Tucker  conditions  characterise  local 
optlna  among  which  is  a  global  optimum.  Modules  are  usociated  with 
partitions  and  a  special  clue,  guillotine  partitions,  are  characterised. 
Branch  and  bound,  partial  enumeration  and  heuristic  procedures  for  finding 
optimum  or  good  guillotine  partitions  are  discussed  and  illustrated  with 
examples. 


1.  INTRODUCTION 


The  one-module  modular  design  problem  was  first  stated  by  David 
Evans  in  [3].  In  the  problem  parts  are  to  be  grouped  into  a  single  module, 
several  of  which  can  then  be  used  to  satisfy  (or  over  satisfy)  the  require¬ 
ments  for  parts  In  a  given  application.  The  objective  is  to  minimise  the 
total  cost  of  the  parts  used  for  all  applications.  The  formal  statement  of 
the  problem  is: 

Minimise  E  c.x.  £  d.y. 
x,y  id  J«J  3  3 

Subject  to 


where 


xi,yj,ci,dj,rij  -  0  J 

and  y^  are  integers 


for  i«I,  jeJ 


I  *  { 1 , 2 , . . .  ,m} 

J  ■  {l,2,...,n} 

Cj  ■  cost  of  part  i 

dj  ■  demand  for  application  j  (an  integer) 

r1J  •  number  of  part  i  units  required  in  application  J  (am  Integer' 
xi  » 'the  number  of  part  1  units  on  the  module  (a  decision  variable) 
yj  -  the  number  of  modules  needed  for  the  Jth  application 
(a  decision  variable) 

The  continuous  modular  design  problem  is  obtained  by  dropping  the  Integer 
requirements  on  and  yj.  A  simplex-like  solutlom .procedure  for  this 

problem  was  c,'"mr'  Shaftel  and  Thompson  in  [t J .  Several  aqpmsdb  procedures  ion 
the  problem  were  previously  given  by  Evans  31,  Charnes  and  Kirby  Cl] ,  and,  more 
recently,  Passey  [6]. 


Evans  in  [4]  ana  Rutenberg  and  Shaftel  In  [7]  have  extended  the 
modular  design  (MD)  problem  to  the  case  where  more  than  one  module  can  be 
used.  This  new  problem  is  called  the  multiple  modular  design  (MMD)  problem. 
Rutenberg  and  Shaftel  formulated  the  problem  for  more  chan  one  market  with 
certain  other  costs  and  devised  a  heuristic  aearcn  procedure  for  locating  good 
integer  solutions.  Recently  Silverman  [10,  11]  has  presented  a  search  proce¬ 
dure  for  solving  the  MMD  problem,  as  formulated  by  Evans,  given  a  solution  to 
the  MD  problem. 

In  the  present  paper  we  shall  extend  the  results  of  our  earlier  paper  [9] 
to  solve  the  continuous  MMD  problem.  It  is  felt  that  solutions  to  the 
continuous  problem  will  lead  to  more  efficient  solutions  of  the  integer 
problem.  We  first  sharpen  the  definition  of  the  MMD  problem  by  relating  it 
to  partitions  of  the  MD  problem.  We  then  define  guillotine  partitions  and 
study  them  in  detail.  Finally  we  develop  branch  and  bound  and  heuristic 
techniques  for  finding  optimal  and  good  guillotine  partitions. 

2.  THE  PRIMAL  AND  DUAL  MMD  PROBLEMS 

The  multiple  modular  design  (MMD)  problem  for  a  single  market  requires 
that  parts  be  grouped  into  several  modules  and  various  combinations  of  these 
modules  be  used  in  applications.  There  is  a  fixed  cost  of  producing  each 
additional  module  and  certain  other  costs  (such  as  handling  costs)  which  are 
a  linear  function  of  the  number  of  modules  needed  in  each  application. 

The  MMD  has  been  formally  stated  [7]  as  follows: 

Minimize  t  [E  c,  x.<k)  I  0,  +  E  F(k)  tilt,  l5k)  *  8 

x.y.P  k  i  1  -*  s  1  5  k  k  j  J 
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Subject  to 


E  Xj<k)  >  Xij  for  leI*  J«J 


t(k)  v<k> 

k  * 


*  ci»  dj»  In*  fe j  £  0  for  ltl»  kfK 


*i 


(k)  ..(k) 


LJ 


J*  *ij*  *j 
integers 


where 


I  *  Cl, 2, ... ,m} 

J  “  Cl, 2 . n} 

K  -  {l ,2 , . . .  ,p} 

end  p  (■  the  number  of  modules)  is  a  decision  variable; 
c^  »  cost  of  r>«?t  i 
d.  *  demand  for  application  j 


■J 
*1J 


■  number  of  units  of  pert  i  required  in  application  J 


bj  ■  cost  of  placing  a  module  in  application  j 
fk) 

F  ■  fixed  cost  of  producing  the  kth  module 
^  «  number  of  units  of  part  i  on  module  k  (a  decision  variable) 


*i 


,(k) 


number  of  units  of  module  k  needed  in  application  J 
(a  decision  variable) 

In  this  paper  we  shall  concentrate  on  the  continuous  version  of  the  >MD 

problem  and,  therefore  will  omit  the  integer  restrictions  on  the  x  and  y 

variables.  For  the  continuous  problem  we  can  make  the  substitutions 

(k)  (k)  (k)  .  (k) 

x'  ■  c.  X'  ,  y;  -  dj 


'i 


’i  *L 


J 


ij 


ci  dj  Xij*  and  bJ  "Vdj 


in  order  to  simplify  the  statement  of  the  problem.  Initially  we  will  restrict 
our  attention  to  the  problem  where  the  number,  p,  of  modules  has  been  fixed. 
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In  this  c«««  if*  c*n  drop  th*  fixed  charge  tent  in  the  objective  function 

since  it  1*  a  constant.  The  problem  then  becomes: 

Minimise  £  t£  x(k)  £  y 00 ]  +  £  £  b  y.(k)  -  g  (IA) 

x,y  k  i  1  j  J  k  J  J  J  p 


subject  to 


r  «<k>  v<k>  >  r 

J  *i  yJ  -  U 

rij*  bj*° 


*(k) 

i 


,(k) 

'i  1 


(2A) 

(3A) 


Later  we  will  show  how  to  calculate  an  upper  bound  on  p  and  give  a 
search  procedure  for  minimising  gp  over  all  possible  values  of  p.  To  do 
this  we  will  have  to  bring  the  fixed  charge  terms  back  into  the  objective 
function. 

We  shall  now  derive  a  simplified  primal  problem,  and  then  use  the  Kuhn* 
Tucker  conditions  to  derive  a  dual  problem  similar  to  the  one  we  found  for  the 
single  module  case  in  [9] . 

THEOREM  1.  In  searching  for  an  optimal  solution  to  problem  (A)  the 

linear  term  in  (1A)  Can  be  Ignored. 

PROOF.  Let  x^k\  for  k«K  be  feasible  for  problem  (A).  If 

(k) 

is  any  number  >  0  then  x  /t^,  for  kcK  is  also  feasible.  By  setting 

•  t  <  1  for  all  k  the  second  term  of  objective  function  for  this  solution 

i«  reduced  while  the  first  term  remains  constant.  By  letting  t  -*  0  the 

linear  term  cen  be  made  arbitrarily  smell,  so  that  problem  (A)  has  no  minimising 
(k) 

solution  with  y  >0.  However  the  objective  function  is  nonnegative  and 
hence  has  an  lnfimum  M  over  the  constraint  set.  If  we  solve  problem  (A)  without 
the  linear  term  we  can  then  use  the  above  transformation  to  obtain  a  solution 
with  value  arbitrarily  close  to  M.  Thus  we  are  justified  in  solving  the  con¬ 
tinuous  problem  without  the  linear  term. 


It  should  ba  remarked  that  for  the  Integer  problem  there  are  additional 

/M 

constraints  of  the  for  a  yjj 1  >  1  eo  that  problem  (A)  has  r  minimum  solution 

that  is  attained  within  the  constraint  set.  In  the  integer  case  the  linear 

tana  cannot  be  ignored,  and  Theorem  1  does  not  hold  for  the  integer  problem. 

In  order  to  solve  the  continuous  MMD  problem  we  can  (as  in  [3,9])  isolate 

(k) 

one  of  the  set  of  optimum  solutions  by  requiring  that  I  y,  ■*  1  for  each  k. 

J«J  J 

fk)  ik) 

But  we  shall  be  even  more  explicit  concerning  the  and  yk  '  variables,  as 
follows. 

Suppose  that  in  the  kth  module  we  choose  to  make  some  of  the  x^  or 
(k) 

yj  variables  equal  to  zero.  In  order  to  derive  a  dual  problem  that  reflects 

this  we  wish  to  record  these  zero  constraints  explicitly.  Define  p  subsets 

I. and  J.,...,J  such  that  I.  c  i  j.  cj  for  keK  and 
1  p  X  p  k  x. 

I  x  J  •  (Ij  X  Jj)  U  (I2  X  J2)  U...U(lp  X  Jp) 

Notice  that  we  do  not  (as  yet)  require  (Ih  x  Jh)  x  (Ifc  x  Jk)  ■  0  for  h  4  k,  i.e., 
we  do  not  require  that  the  sets  1^  x  Jk  partition  I  x  J.  We  now  call  the 


Subject  to 


Minimize  £ 

X  k*K 

,■  *(k)  -  . 
tc\  1  ^ 

(IP) 

I  »<k)  y.(k>  ■ 

keK  1  J 

Zij  "  rij’  i,X»  j*J 

<2P) 

£  y^  •  1 

3 

keK 

(3P) 

•  o 

keK,  iel-lk 

(4P) 

yjk)  .  0 

keK,  JeJ-Jk 

(5P) 

z  x(k)  v(k) 

ij’  *i  *  yj 

>  0  keK,  iel,  JeJ 

(6P) 

where  z^  is  e  surplus  variable. 
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THEOREH  2.  An  optimal  solution  to  problem  (?)  exists. 

PRC')?.  Problem  (P)  is  to  minimite  a  continuous  function  over  a  compact 
sat  and  by  a  well-known  theorem  of  mathematics  has  an  optimal  solution. 

In  this  paper  we  shall  indicate  some  algorithms  and  some  heuristic 
procedures  for  finding  exact  and  approximate  solutions  to  the  continuous  MMD 
problem. 

LEMMA  1.  If  there  are  no  zero  rows  or  zero  columns  in  the  requirements 
matrix,  R,  then  in  any  feasible  solution,  for  each  i  there  is  at  least  one 
k  with  x<k>  >  0,  and  for  each  J  there  is  at  leaBt  one  k  with  yjk^  >  0. 

PROOF.  Since  2  x,^  y|k^  >  r .  .  >  0  for  a  feasible  solution,  it 
1-  *  J  “  J 

*  He}  (It} 

follows  that  there  is  at  least  one  k  with  x^  '  y^  >  0. 

fit) 

Let  be  variables  associated  with  the  constraints  (2P)  and  let  g 

be  variebles  associated  with  constrrints  (3P).  The  Kuhn-Tucker  conditions 
associated  with  the  primal  problem  can  be  shown  to  be: 


E 

^k 

x,j 

v(k)  -  1 

yJ 

for 

keK, 

UIk 

(Id) 

E 

U1k 

x.) 

(k)  (k) 

x^  -  gv 

for 

kcK, 

J"k 

(2B) 

XtJ  -0 

for 

HI, 

jeJ 

(3B) 

xl)i° 

for 

id, 

jeJ 

(4B) 

LEMMA  2.  For  any  pair  of  dual  feasible  ‘olutions,  2  xfk^  ■  g^ 

leIk 

for  kcK.  K 

(k) 

PROOF.  Multiply  equation  (IB)  by  x^  and  sum  over  iel^;  obtaining 

E  E  yjk)  x,(k)  -  E  x.(k)  for  keK. 

l>  >  uik 

Similarly,  multiply  (2B)  by  y^  and  sum  over  jeJ^  and  use  (3P)  to  obtain 
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-  ■ - **■ 


z  z 
J«Jk  i«Ik 


X,,  y 


(k) 


'ij  'J 


z 


g(k)  for  k«K 


Combining  these  two  equations  gives  the  desired  result. 

(k) 

LEMMA  3.  For  any  pair  of  dual  feasible  solutions  Z  g  *  g  . 

kcK  p 

PROOF.  This  follows  from  Lemma  2  and  (IP). 

By  analogy  to  the  continuous  MO  problem  (see  [9])  and  classical  linear 
programming ,  we  shall  create  a  dual  problem  from  constraints  (IB),  (2B),  and  (4B) 
and  the  objective  function.  The  dual  problem  is: 


Subject  to 


Maximize  Z  Z  X.  .  r 

X  i  1 

*ij  J 


Z 

J«Jk 


(k) 
ij  'J 


X,,  y 


1 


for  ktK,  iel^ 
for  k«K,  jeJ^ 


(ID) 


(20) 

(30) 


Xtj  >  0  for  UI,  jcJ  (4D) 

Constraints  (3B)  can  be  interpreted  as  complementary  slackness  conditions. 
They  will  be  forced  to  hold  by  the  solution  procedures  wc  shall  present. 

Note  that  for  each  module,  the  constraints  for  fixed  k  have  a  similar 
format  to  that  of  the  dual  problem  of  the  continuous  MD  problem  (see  19] )  but 
restricted  to  the  1^  x  area  of  the  matrix  R.  Many  properties  of  the  con¬ 
tinuous  MD  problem  will  carry  over  to  the  continuous  MMD  problem. 

LEMMA  4.  For  any  pair  of  feasible  primal-dual  solutions  (whether  non- 
negative  or  not)  we  have 


g 


P 


Z  Z  X 
i  ) 


ij  zij 


sum  over  all  1  and  j  and  use  (ID)  to  show: 


PROOF.  Multiply  <2P)  by  X^, 

im,.  x.<k>  yjk)  - 
ljk  1  J  1  J 


■ij 


r  z 
i  J 


ij  ij 


£P  ‘ 


As  in  the  proof  of  Lemma  2  we  also  have 


EIEL. 

ijk 


x(k>  v<k> 
i  yj 


Z  Z  x 
k  i 


(R) 

i 


8 


P 


which  proves  the  lsmma. 

LEMMA  S.  (Complementary  slackness).  For  any  pair  of  dual  feasible 
solutions  gp  •  f  if  and  only  if  X^  z^  *  0  for  iel,  )«J. 

PROOF.  The  proof  follows  directly  from  Lenma  A  and  the  fact  that  both 
Xj^  and  z^  are  nonnegative. 

THEOREM  3  (Duality  Theorem).  A  necessary  condition  that  x^k\  yjk^  be 
optimal  for  problem  (?)  is  that  there  exists  a  solution  to  problem  (D) 

end  that  g  ■  f  . 

P  P 

PROOF.  Using  Letama  l,  we  can  show  that  the  Arrow-Hurwicz-Uzawa  constraint 
qualification  [3,  p.  102]  holds  for  problem  (P).  Therefore  at  the  optimum 
solution  to  problem  (P)  there  exists  a  solution  to  the  Kuhn-Tucker  problem, 
since  problem  (D)  together  with  X^  z^  ■  0  (which  is  true  if  and  only  if 
gp  ■  f  from  Lasasa  3)  make  up  the  Kuhn-Tucker  conditions.  These  conditions 
swat  hold  at  the  optimum  because  of  the  constraint  qualification  [5,  pp.  103-106] . 

The  MD  problem  had  a  unique  solution  so  the  duality  theorem  for  that  case 
was  stronger  then  here.  In  the  MMD  problem  there  are  many  pairs  of  solutions 
satisfying  the  conditions  of  Theorem  3,  and  these  correspond  to  IqcjI  ogtj^a 
of  g^.  Among  these,  of  '-curse,  is  a  global  optimum. 


3.  CHARACTERIZATION  OF  LOCAL  OPTIMA 

In  the  present  section  we  shall  give  some  properties  of  local  optima 
for  tin  case  of  a  fixed  number,  p,  of  modules.  It  will  be  shown  that  a  local 
optimum  to  the  problem  can  be  found  by  solving  p  MD  problems. 
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DEFINITION.  By  «  partition  of  a  MKD  problem  with  data  matrix  R  (bare* 

(k) 

aftar  called  problem  R)  into  p  problem!  with  data  matrices  R'  (hereafter 

called  problem  R^)  for  k«R  we  mean  R  ■  I  R^\  i.e. 

kcK 

r. .  »  £  rfk^  where  r^  >  0. 

kcK  J 


An  integral  partition  is  one  that  also  satisfies, 

rj^*  ■  ct  dj  j^k*  where  >  0  is  an  integer 


We  shall  regard  each  R 


(k) 


as  defining  a  MD  problem  with  variables  x. 


(k)  ..(*»;) 


and  X 


(k) 

IJ  ‘ 

THEOREM  (A)  Given  an  optimal  primal  dual  solution  x^k\  y^  and 

(k) 

X,  ^  to  problem  R  a  partition  R  ■  £  R  can  be  defined  in  such  a  way  that 
x^\  y^  and  X^  ate  optimal  primal-dual  solutions  to  the  MD  problem  R^k\ 


(B)  Given  a  partition  R  *  £  R 


primal  solutions  for  the  MD  problem  R 

,<k) 


(k) 

1 

<k> 


let  x<k)  and  y^  be  optlsial 

Then,  if  it  is  possible  to  choose 
(k) 

X^  for  all  keK  as  optimal  dual  solutions  to  R  for  all  kcK,  then 


and  y 


(k) 

J 


give  a  local  optimum  for  problem  R. 


(C)  If  g^  is  the  value  of  problem  R^  and  the  condition  of  (B) 


holds  then 


«p  •  2  8 

p  keK 


(k) 


PROOF.  (A)  From  (2P)  we  have 


r  ,<k>  v<k>  >  r 

r  i  yJ  - 


.  .  (k)  .  (k)  (k) 

hence  we  can  choose  r^7  <  x^  7  yj  so  that 


P 

£  ri « 

k-1  1J 


(k) 


r^  for  each  1  and  j 


in  at  least  one,  and  perhaps  many  different  ways.  If  z^  *  0  then 


t. <»,<«- f  tl 


so  .<$>  «  0  for  problem  R^. 


Hence  X^ j  >  0 


.<k> 


implies  s . ,  »  0  and  complementary  slackness  holds.  The  primal  and  dual 
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(k) 

conditions  for  problem  R  are  as  follows: 


(1) 

,00  »00  >  r00 
*1  yj  ^  rij 

for  ielk,  JcJ, 

(2) 

£  y(*0  m  i 

j*k  J 

for  kcK 

(3) 

I  x<k>  -  g<k> 
itlfc 

for  kcK 

(4) 

r  \  v<k>  •  1 

jik  1J  J 

for  kcK,  iclk 

(5) 

E  X4.  x<k)  -  g(k) 

i«Ik  J 

for  kcK,  JcJ 

(6) 

£  £  X  r<k)  -  f 

i«Ik  Wk  iJ 

(k)  for  kcK 

(7) 

(k)  (k)  .  0 

*i  .  yj  .  xtj  >  0 

for  keK,  iclk 

Uo  have  already  shown  that  (1)  holds.  Equation  (2)  is  the  same  as  (3P);  (3) 

and  (6)  are  true  by  definition;  (4)  Is  (2D);  (5)  is  (3D);  and  (7)  Is  included 

(k)  (k) 

In  (4P)  and  (4D).  Hence  ,  yj  and  j  are  primal-dual  solutions  to 
problem  R^  . 

(B)  If  for  a  partition  £  ■  r  ve  can  choose  ■  X^ j  for  ail 

i  and  J,  then  the  steps  in  the  proof  of  (A)  are  reversible.  We  omit  the  rest  of 
the  details. 


(C)  If  we  sum  (6)  over  k 

and  the  partition  formula  we  have 

£  f^  »  £  £  £ 

k  kcK  Ulk  JcJ 


k 


and  use  the  duality  theorem  for  problem 


£ 

kcK 


f 

P 


R<k> 


completing  the  proof  of  the  theorem. 

By  means  of  this  theorem  we  see  that  one  way  of  finding  local  opt lata  to 

(k) 

the  HKD  problem  is  to  partition  R  into  HD  problems  R'  and  then  use  the 
algorithm  of  (9]  to  solve  R^. 


If  it  happens  that  the  seme  dual  variables 
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00 

can  be  used  for  each  problem  R  Chen  we  have  found  a  local  optimum  to  Che 
NMD  problem.  If  we  can  enumerate  all  local  optima,  we  can  choose  Che  smallesc 
as  Che  global  opClmum.  There  are,  of  course,  an  infinite  number  of  arbicrary 
particions.  We  shall  concenCraCe  on  incegral  partitions,  since  there  are 
only  a  finite  number  of  them. 

DEFINITION.  By  a  guillotine  partition  with  p  pieces  of  R  we  shall 

mean  a  partition  of  I  x  J  into  p  "rectangular"  subsets 

I  x  J  -  (I.  x  J.)  U...U  (1  x  J  ) 

11  P  P 

where 

(1^  x  J^)  0  (I^  x  J^)  *  0  for  h  +  k 
(k ) 

and  a  partition  R  ■  E  R  where 

keK 

rij)  “  rij  if  (1,J)  G  lk  *  Jk 

-  0  if  (i,j)  «  Ih  x  J(i  and  h  4  k. 

By  a  guillotine  row  partition  we  mean  a  guillotine  partition  with 
*  J  for  keK.  Similarly  by  a  guillotine  column  partition  we  mean  a  guillotine 
partition  with  1^  »  I  for  keK. 

Clearly  guillotine  partitions  are  integral  and  there  are  only  a  finite 
number  of  them.  Also  any  guillotine  partition  can  be  constructed  stepwise  by 
row  and-or  column  guillotine  partitions  applied  alternately  to  a  series  of  sub¬ 
problems  . 

(k) 

THEOREM  5.  Suppose  problem  R  is  divided  into  p  rectangular  pieces  R 
by  a  guillotine  partition.  If  and  xj^  are  primal-dual  optimal 

solutions  for  R^  then  ,  y^  and  X  *  £  > .  ^  give  a  local  optimum 


to  problem  R. 
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(k)  (k) 

PROOF.  Suppose  R  consists  of  the  matrix  defined  by 


.00 

ij 


■ij 


r(k)  -  0 

ij 


for  (i ,  j )  €  Ik  x  Jk 
for  (i,j)  i  Ifc  x 


,(k) 


Then  if 


where  Ik  c  I  and  Jk  c  j  define  the  rectangular  piece  R 
id^  we  have  x^  »  0  and  if  j<Jk  ve  have  ■  0.  Hence  x^  ■  0  if 

(i,J)  i  Ik  x  Therefore 


(k) 

lU 


0  if  (k,J)  i  Ik  x  Jk»  It  follows  that  if  we 


define  X 


ij 


-  \ 


(k) 


it  will  satisfy  the  dual  conditions  (2B)-(4B).  Hence 


k-1 


(k)  (k) 

by  Theorem  4  x^  and  y  give  a  local  optimum  for  problem  R. 


LEMMA  6.  Let  Ik  x  Jk  be  a  rectangular  piece  of  a  guillotine  partition 

al  vi 

"ij* 


(k)  (k) 

and  let  g  be  the  optimal  value  of  R  considered  as  an  MD  problem. 
Then  g(k)  >  £ 


(i,J)*IkxJk 


PROOF.  Since  X^  »  1  for  (i,j)  e  Ik  x  Jk  is  always  dual-feasible  and 
(k)  (k) 

g  ■  f  by  the  duality  theorem,  the  result  is  obvious. 


.00 


DEFINITION.  A  sub  problem  R 

£  r  i  1  * 

(i,J)cIkxJk  lJ 


(k) 


defined  on  Ik  x  Jk  is  tight  if 


THEOREM  6.  Consider  the  MMD  problem  with  all  fixed  costs  equal  to  •>, 

(a)  For  guillotine  partition  gp  >  gp+^» 

(b)  Unless  all  pieces  of  a  guillotine  partition  with  p  pieces  are  tight, 
there  is  an  additional  row  or  column  guillotine  partition  that  can  be  made  such 

«>«  gp  >  8p+1. 

PROOF,  (a)  The  guillotine  partition  with  p+1  pieces  is  obtained  from 
that  with  p  pieces  by  applying  a  row  or  column  guillotine  partition  to  one 
of  the  original  pieces.  Suppose  Ij  x  is  to  be  divided  into 
<lux  Jj)  U  (1^2  *  Jj)  by  a  row  partition  (the  proof  for  column  partition 


0 
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is  similar).  Then  fche  y  vector  that  was  optimal  for  1^  x  together  with 
the  1^^  components  of  x  will  be  feasible  for  1^  x  J^.  Similarly  the  sane  y 
vector  and  the  1^  components  of  x  will  be  feasible  for  1^  x  Jj.  Hence 
(a)  follows. 

(b)  If  the  1^  x  piece  for  the  guillotine  partition  with  p  pieces 
is  not  tight  then  some  cell  (i,J)  c  1^  x  has  >  0.  By  making  the  next 
partition  isolate  either  the  row  or  the  column  containing  (i,j)  we  can  assure 

eh*t  8p>8p+i- 

We  give  next  an  example  of  a  problem  for  which  the  optimum  partition  is  not 
a  guillotine  partition.  Consider  the  problem 


and  suppose  we  require  p  ■  2.  The  sum  of  the  elements  in  R  is  27,  and  the 
reader  may  easily  verify  that  every  guillotine  partition  has  g2  >  27,  However 
the  following,  non  guillotine,  partition  achieves  g2  *  27: 


1 

2 

3 

i° 

l 

2 

T 

— 

1 

l 

R  « 

2 

3 

A 

-  |o 

l 

2 

+ 

2 

2 

2 

3 

A 

5 

!o 

l 

,._2 

3 

3 

3 

-r<1>+r<2) 

The  optimum  solution  to  is 

x  -  (3,  3,  3),  y  -  (0,  |), 

(21 

and  the  optimum  solution  to  R  'is 

x  ■  (3,  6,  9),  y  •  (-|,  -j,  |) 

One  optimum  dual  solution  is  *  1  for  all  i  and  j,  and  there  are  others. 
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IC  is  interesting  to  note  that  the  vertical  guillotine  partition  that 
defines  coluan  1  as  one  subproblem  and  columns  2  and  3  as  the  other  sub- 
problem  has  a  value  g^  ■  27.58  which  is  only  slightly  more  than  the  optimum 
of  27.  We  have  generally  found  that  guillotine  partitions  will  give  local 
optimum  values  close  to  the  optimum  when  they  do  not  themselves  provide  the 
optimum. 

in  the  rest  of  this  paper  we  shall  restrict  ourselves  to  guillotine 
partitions.  Other  reasons  why  they  are  advantageous  is  that  row  guillotine 
partitions  restrict  the  number  of  different  parts  that  appear  on  a  given 
module,  while  column  guillotine  partitions  restrict  the  number  of  different 
modules  that  appear  in  a  given  application.  Both  of  these  conditions  are 
probably  desirable  from  a  manufacturing  point  of  view,  even  though  they  are 
not  specifically  included  in  the  objective  function. 

We  intend  to  take  up  the  question  of  non-guillotine  partitions  in 
another  paper. 

4.  SOLUTION  PROCEDURES  FOR  GUILLOTINE  PARTITIONS 

Searching  successive  MD  problems  of  a  MMD  problem  is  made  easier  by 
the  ability  to  use  a  type  of  parametric  programing.  Given  the  optimal 
solution  to  a  MD  problem,  we  can  readily  determine  a  ge jd  initial  feasible 
solution  for  any  subproblem  formed  via  a  guillotine  partition,  as  follows. 

Let  x^,  id,  y^,  J«J  be  the  optimal  solution  to  a  particular  mxn  MD  problem. 
Now  assume  that  we  partition  the  problem  into  two  subproblims  made  of  colunns 
of  the  original  problem  (the  same  reasoning  would  apply  for  rows).  Let 

■  {the  set  of  columns  which  are  posirive  in  moduli,  k},  k  ■  1  or  2. 


. 


*  ££££•-.::  jaar-imsMag.  W 


'  '****•  ,  «  **  ?i  .msh^tif.  C«-*  =**«#>%  VWWP*'**^  g‘Og*,W*^?n» *»«W ' 
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tod 


yjk)  ■  yj /  £  yj  j«Jk  «nd  k  -  1  or  2 

J«Jk 


.<k) 


(k) 


.<k) 


Let  Che  cell*  which  determine  be  eight  end  by  increasing  the  right 

hand  sides  of  other  constraints  force  these  cells  to  enter  the  basis  until 
a  tree  basic  solution  is  formed  with  r^*  ■  x ^  yjk^  for  (i,J)  in  the  tree. 

THEOREM  7.  (A)  The  solution  described  above  is  primal  feasible. 

(B)  If  m  and  n  are  >  1  and  the  original  problem  is  nondegenerate 
the  sura  of  the  objective  functions  for  the  initial  primal  feasible  solutions 
to  the  subproblems  is  less  than  the  optimal  solution  for  the  original  problem. 


PROOF.  (A)  Let  rlT/y<k)  -  max  r^/yj10 

J*Jk 


than,  xjk)  yjk)  -  rlT/y*k)  •  yjk)  >  rlj/yjk)-yjk)  -  r^  for  all  i,j  and  k. 

(B)  Let  g  be  the  optimal  solution  to  the  original  problem.  And  h^ 

and  hg  be  tba  objective  functions  associated  with  the  initial  primal  feasible 

solutions  of  the  two  subproblems.  Then. 

g  •  I  E  x.  y..  We  also  know  that  x.  ■  max  r.  /y  for  all  i, 
UI  jeJ  J  jeJ  J  J 

so  that  g  ■  E  E  y.  •  max  r  /y  . 

i«I  jeJ  J  JeJ  J  J 


And, 


h1  + 


£  y  oax  r../y.  +  E  y.  max 
j«Jl  J  jeJt  J  jd2  J  jcJ2 


ru/yi 


But  for  m  and  n  >  1  (and  a  nondegenerate  problem)  there  Is  at  least  one 

row  in  R  which  contains  only  one  tight  cell  so  that 

oax  r.  /y.  >  max  r../y.  . 
j«J  J  j*Jk  J 

Hence,  g  >  hj  +  hj. 


A  solution  for  a  subproblem  formed  in  the  above  manner  is  usually 
a  very  good  Initial  primal  feasible  solution.  It  Is,  of  course,  an  upper 
bound  on  the  value  of  the  objective  function  for  that  suhprobiem.  A  lover 


bound  to  the  objective  function  value  of  a  subproblem  is  £ 


ieIk  JCJk 


.(*) 

ij 


It  is  important  to  note  that  the  upper  bound  calculation  is  the  first  (and, 
probably,  moat  time  consuming)  step  of  calculating  the  exact  optimum  for  a 
*ubproblem.  Because  of  this,  it  seems  likely  that  for  any  branch  and  bound 
emaaeration  procedure  it  will  pay  to  calculate  exact  solutions  rather  than  use 
approximate  lower  bounds.  It  is  possible  to  enumerate,  explicitly  or  implicitly, 
all  possible  ways  of  partitioning  the  problem  into  p  modules  (bounds  on  the 
value  of  p  will  be  discussed  in  the  next  section).  The  organization  of  such 
an  enumeration,  not  the  calculation  of  lower  bounds,  is  what  would  make  such 
a  complete  enumeration  difficult.  Because  of  these  difficulties  we  will  present 
some  properties  which  sould  make  useful  tools  for  obtaining  good  heuristic 
solutions. 

Given  a  solution  to  a  particular  MD  problem,  the  greatest  value  by  which  we 

can  reduce  the  value  of  the  objective  function  is  £  £  z  If  we  form  a  partition 

i  j 

with  a  subproblsm  with  a  single  positive  column,  then  the  objective  function 
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will  reduce  by  at  least  £  z.  ,  for  that  row.  Since  m  is  generally  larger 

J  3 

than  n,  vertical  partitions  seem  the  most  appropriate.  Another  rationale  for 
using  only  vertical  partitions  can  be  made  by  counting  the  number  of  new  tight 
cells  which  would  be  caused  by  vertical  versus  horizontal  partitioning.  Once 
a  vertical  partition  is  made,  subdividing  the  new  subproblems  into  more  vertical 
partitions  is  even  more  appropriate.  If  a  total  of  n  partitions  were  to  be 
made,  n  modules  each  satisfy  the  needs  of  exactly  one  application  would  yield 


the  optimum  value  of  £  I  r 

1  J 


U* 


Given  the  value  of  p,  or  an  estimate  for  that  value,  and  using  vertical 
partitioning  only,  it  becomes  necessary  to  intelligently  group  applications  to 
be  satisfied  by  certain  modules.  (Note  that  the  techniques  discussed  here  will 
also  apply  to  horicontal  partitioning.)  This  fact  would  mean  that  we  would 
want  to  choose  applications  to  group  whose  parts  requirements  have  approximately 
the  same  ratios.  We  would  also  like  to  group  applications,  with  requirements 
which  are  difficult  to  fulfill,  into  the  smallest  groups.  Three  possible 
measures  which  may  be  used  are: 

1).  The  variance  (or  absolute  difference)  along  each  column  is  an 
Important  measure.  Since  we  are  dealing  with  ratios,  we  would  want  a 
measure  of  the  variance  of  the  normalized  row  entries.  The  measure  of  this 
type  that  we  shall  use  is 


£  |rtj/£  rtj 


2).  The  above  value  would  be  useless  without  some  sort  of  indication 


of  the  size  order  of  the  entries  of  a  column.  Ideally,  of  course,  we  would 
compare  each  pair  of  columns  to  see  how  closely  they  fit.  An  easier,  but 
lass  informative,  method  is  to  compare  each  column  with  a  given  standard  order 
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Tha  sum  of  the  absolute  difference  between  the  order  and  the  standard  order 
will  be  used  as  a  measure. 

3).  Clvan  that  we  make  only  vertical  partitions  the  number  of  columns 

In  a  given  module  Is  likely  to  be  on  the  order  of  magnitude  of  n/p.  In  order 

to  determine  which  colisnns  should  be  In  the  smaller  partitions  we  look  at 

£  z  for  the  single  module  solution.  Columns  with  larger  values  of  £  z 
1  J  1 

would  be  in  partitions  which  satisfy  the  least  number  of  applications. 

We  consider  the  effectiveness  of  the  three  heuristics  for  a  specific 

example  in  Section  6. 

5.  DETERMINING  THE  NUMBER  OF  MODULES 

As  mentioned  earlier,  n  la  an  upper  bound  on  the  value  of  p  while  1 
la  a  lower  bound.  The  solutions  for  these  two  cases  are  quite  easy  to  find; 
the  difference  between  the  two  solutions  (excluding  fixed  costs)  Is  exactly 


£  £  i 

i  J 


ij 


where  z^  is  the  surplus  for  the  case  where  p  “  1.  Reduction 


in  the  value  of  Z 


(k) 


as  p  increases  soem3  to  decrease  at  least  linearly. 


Aa  more  experience  in  solving  MMD  problems  is  obtained,  good  estimates  on  p 

should  become  available.  At  present  a  good  estimate  (assuming  is  an 

(o)  Z(1) 

Increasing  function  of  p)  is  the  largest  value  of  p  such  that  Fvp'  <  . 

During  a  search  for  the  exact  value  of  p ,  an  upper  bound  which  is  often  better 

than  n,  as  well  as  a  means  for  eliminating  certain  possible  values  between  ]  and  n 

can  be  achieved.  These  two  procedures  are  outlined  in  the  next  two  theorems. 

THEOREM  8:  Let  *  £  Z  z .  .  when  q  modules  are  produced.  Then  an 

i  j 

upper  bound  on  p  is  the  largest  Integer  such  that: 
z(q)  .  ^  F(k)  >  0 
k*q+l 


*/n**n. 


-19- 

fk) 

where  F  It  the  fixed  cost  of  the  kth  additional  module. 

PROOF.  By  definition  Z^  is  the  maximum  we  can  reduce  the  value  of 
the  objective  function  by  producing  additional  modules.  As  soon  as  the  total 
coat  of  additional  modules  surpasses  this  amount,  adding  a  module  cannot  be 
profitable. 

THEOREM  9:  Given  Z^  and  6,  the  value  of  the  best  minimum  solution 
found  so  far,  calculate  c  as  the  largest  integer  such  that 

J  Jr  +  z*q^  +  I  <  e 
i  J  iJ  k-1 

Than  any  value  of  k  from  (and  including)  c  to  q  cannot  be  the  optimum. 

PROOF.  Since  Z^  is  a  nonincreasing  function  of  q,  it  is  a  lower 
(c) 

bound  on  the  value  of  Z  ,  c  <  q.  This  implies  that  a  lower  bound  on  the 

value  of  the  objective  function  for  any  value  of  c  <  q  is: 

E  E  r .  .  +  Z^  +  E  Fk.  We  need  only  consider  values  of  c  which  have  a 
i  j  k-1 

lower  bound  less  than  the  optimum  found  thus  far.  Hence  the  theorem. 

6.  EXAMPLE 

We  will  give  an  example  problem,  enumerate  all  possible  guillotine 
partitions,  and  Indicate  the  values  of  the  heuristic  measures  defined  in 
section  4,  The  problem  we  will  use  1b  symmetrical  so  that  only  vertical 
partitions  need  be  considered. 


1 

2  1  3 

i 

r. 

4 

J 

9  •  6 

12 

10 

3 

6  '  16 

20 

25 

» — 

1 

- 

»-  •  - 

l 

12  !  20_ 

U2-5ii3J 

4 

10  25 

mm-m  - 

23 

36 

Er.  -  11  39  70  81  98 

J  J 

E  E  r. .  -  299 
l  J 

_ Li 
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For  P  -  1  g<l)  -  g{l)  -  367.36.  Z(1)  -  68.36.  For  p  -  2  there  are 

13  possible  vertical  guillotine  oartitions.  In  the  following  table  1,2, -3,4, 5 
means  subproblam  one  has  r^j  for  columns  1  and  2  and  zeros  in  columns  3,  4, 
aad  3,  while  subproblem  2  has  columns  3,  4,  and  5  at  the  value,  r^,  and  zero's 
In  columns  1  and  2. 

Partition  Sub-module  1  Sub-mo iule  2  Total 


Optimun  Excess  Optimum  Excess  Optimum  Excess 


1  -  2, 3, 4, 5 

11 

0 

343.0 

55 

354.0 

55 

2  -  1,3, 4,5 

39 

0 

299.1 

39.1 

338.1 

39.1 

3  -  1,2, 4, 5 

70 

j 

285.1 

56.1 

355.1 

56.1 

4  -  1,2, 3, 5 

81 

0 

266.9 

48.9 

347.9 

48.9 

5  -  1,2, 3, 4 

98 

0 

247.0 

46.0 

345.0 

46.0 

1,2  -  3,4,5 

57.4 

7.4 

276.6 

27.6 

334.0 

35.0 

1.3  -  2,4,5 

91.1 

10.1 

262.5 

44.5 

353.6 

54.6 

1.4  -  2,3,5 

100.7 

8.7 

245.0 

37.9 

345.7 

46.7 

1,5  -  2,3,4 

118.8 

9.8 

226.0 

36.0 

344.8 

45.8 

2,3  -  1,4,5 

126.9 

17.9 

221.9 

31.9 

348.8 

49.8 

2.4  -  1,3,5 

134.6 

14.6 

198.5 

19.5 

333.1 

34.1 

2,5  -  1,3,4 

160.2 

23.2 

185.4 

23.4 

245.6 

46.6 

3,4  -  1,2,5 

165.6 

14.6 

179.9 

31.9 

345.5 

46.5 

3,5  -  1,2,4 

176.4 

8.4 

153.0 

22.0 

329.4 

30.4 

4,5  -  1,2,3 

200.3 

21.3 

145.4 

25.4 

345.7 

46.7 

For  p  •  2 

then 

the  optinuss  puti! 

•■ion  i s  3,5  -  1 

.2.4  ,<»■ 

329.4, 

Z2  -  30.4. 

From  the  enumeration  of  p  »  2  we  can  enumerate 

all  possible 

solutions  for 

p  •  3  and  p 

•  4. 

For  p  *  3 

the 

best  partitions 

are  1,2  -  3,5-4,  g 

<3)  -  314.4 

E<3>  -  15.4. 

For 

p  »  4  the 

best  partitions  are 

1.2  -  3-4-5, 

,<4)- 

306.4, 

Z<*>  -  7.4, 

Note 

that  the  approximation  for  the 

reduction  in 

zoo 

given  in 

the  last  section 

seems  to  be 

very 

accurate. 

Finally,  we  will  show  the  heuristic  measures  discussed  in  section  4. 


m 


-21- 


for  the  single  module  case: 


Sun 


0 

1 

1.2  ]  4 

2 

1 

0 

6.5  3 

1.2 

6.5 

1.4  ;  .8 

0 

4 

3 

.8  {  0 

7 

:  2 

8 

1 

o  i  r 

0 

8.2 

18.5 

9.9  14.8 

17 

Normalised  value  of 


r 


1J 


Ave, 


.091  ’  .051 

.043 

.012  I  .041 

-  -  J*.  .  .  . 

... 

.148**]  .102  * 

.182  ]  .231 

.086 

i 

.273  *  .154 

.229 

.247  j  .225 

-  _  . 

j 

.091  .308 

.286 

.309  i  .235 

t  ■  •  i . H - 1 


.364  _^_*256  j  .357  .284  j_  .368 

.2  .2  .2  .2  .2 


Sun  of  the  absolute 
difference  fron 

tbe  mean  .473  .528  .543  .480  .515 


Order  (ties  are  given  the  average  value  of  their  possible  positions): 


J 


Calculating  tha  turn  of  th«  absolute  value  of  the  differences  between  pairs  of 
colunns  we  gat 


pair 

value 

pair 

value 

1,2 

7 

3,4 

2 

1,3 

5 

3,5 

2 

1,4 

7 

4,5 

4 

1,5 

3 

2,3 

4 

2,4 

2 

2,5 

6 

Note  that  the  opt laun  solution  for  p  •  2  is  3,5  -  1,2,4  where  two 
rows  which  are  similar  along  the  order  and  variance  measures  (3  and  5)  with 
large  t^'a  have  been  grouped  Into  a  smaller  module.  Although  these 
measures  are  not  perfect,  they  do  give  an  a  priori  indication  of  which 
partitions  are  reasonable. 

As  a  second  example,  we  will  solve  tho  problem  presented  by  Evans 
in  [3] ,  for  p  •  2. 


15 

23 

44 

13 

13 

0 

IS 

17 

35 

34 

12 

22 

77 

65 

101 

For  this  problem  Evans  in  [4]  gives  a  solution  of  270.  Silverman  [11]  gives 
an  optimal  solution  of  269.06  via  a  search  technique.  (Beth  solutions  are 
for  non-guillotine  partitions.)  He  shall  show  one  possible  vertical  guillotine 
portion,  chosen  from  the  magnitude  of  E  t. .  from  the  single  module  solution; 


\ 

i 


*?£**«*:wv  •***,., 


^WSl 
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3 

51.0 

In  chit  cate  g^  -  166.9,  g*  *  101, 
few  minutes  by  hand  via  the 
subproblems. 

7.  CONCLUSIONS 

In  this  paper  wa  have  extended  the  MD  solution  of  [9]  to  the  tMD 
problem.  We  have  characterised  the  local  optimum  and  identified  a  useful 
subset  of  these  optimum  aa  elements  of  a,  possibly,  exhaustive  search.  In 
solving  the  MID  problem  by  the  methods  outlined  in  this  paper  several  advantages 
axe  available,  as  follows. 

(1) .  The  solution  is  based  on  solving  the  MD  problem  which  can  be 
solved  quickly.  Use  of  previous  solutions  for  each  MD  subproblem  will  start 
the  Shaftel-Thompson  MD  algorithm  at  an  initially  good  primal  feasible  solution. 

(2) .  We  have  stated  conditions  that  limit  the  search  for  p.  Also,  any 

branch  and  bound  procedure  for  searching  for  an  optimum  solution  to  a  given  p, 
generates  a  tree  which  can  be  used  during  the  solutions  procedure  for  a  new 

value  p '.  This  means  that  rather  than  solving  several  problems  from  the 
beginning  wa  may  merely  extend  the  tree  for  each  new  value  of  p. 

(3) .  The  solutions  for  p  ■  1  and  p  ■  n  can  be  found  quickly  and 
give  a  good  initial  upper  bound  on  the  value  of  the  objective  function. 

(4) .  There  exist  many  good  local  solutions  to  the  problem  so  that  a 


these  values  are: 

column  1  2 

£.tJ  30.3  14.5 

We  shall  solve  the  partitions  1,  2  -  3. 

2 

g  •  267.9.  This  solution  was  found  in  a 
Shaftel-Thompson  algorithm  applied  to  the 


heuristic  search  will  tend  to  find  a  reasonable  solution. 
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